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ABSTRACT 

« ; 

. ^ 1 Context. Starting from accurate magnetohydrodynamic flux rope equilibria containing prominence condensations, we 

initiate a systematic survey of their linear eigenoscillations. This paves the way for more detailed prominence seismology, 
which thus far has made dramatic simplifications about the prevailing magnetic field topologies. 

Aims. To quantify the full spectrum of linear MHD eigenmodes, we require knowledge of all fiux-surface localized 
modes, charting out the continuous parts of the MHD spectrum. We combine analytical and numerical findings for the 

' , ' continuous spectrum for realistic prominence configurations, where a cool prominence is embedded in a hotter cavity, 

Ph I °r where the flux rope contains multiple condensations supported against gravity. 

■ Methods. The equations governing all eigenmodes for translationally symmetric, gravitating equilibria containing an 
axial shear flow, are analyzed, along with their flux-surface localized limit. The analysis is valid for general 2.5D 
equilibria, where either density, entropy, or temperature vary from one flux surface to another. We analyze the intricate 
mode couplings caused by the poloidal variation in the flux rope equilibria, by performing a small gravity parameter 

Q . expansion. We contrast the analytical results with continuous spectra obtained numerically. 

$_i ■ Results. For equilibria where the density is a flux function, we show that continuum modes can be overstable, and we 

' present the stability criterion for these convective continuum instabilities. Furthermore, for all equilibria, a four-mode 

coupling scheme between an Alfvenic mode of poloidal mode number m and three neighboring (m — 1, m, m + 1) slow 
i— i, modes is identified, occurring in the vicinity of rational flux surfaces. For realistically structured prominence equilibria, 

this coupling is shown to play an important role, from weak to stronger gravity parameter g values. The analytic 
predictions for small g are compared with numerical spectra, and progressive deviations for larger g are identified. 
Conclusions. The unstable continuum modes could be relevant for short-lived prominence configurations. The gaps 
, created by poloidal mode coupling in the continuous spectrum need further analysis, as they form preferred frequency 

ranges for global eigenoscillations. 

Q\ • 

■ Key words, solar physics, solar prominences - Instabilities - Magnetohydrodynamics (MHD) - Plasmas 

2 ' 1- Prominence seismology from MHD modeling 

Motivate d by the recent interest an d progress made in diagnosing oscillations in solar prominences, we have started 
an effort (jBlokland fc Keppensll20i"ll hereafter Paper I) to close the gap between observed prominence properties, and 
the models used to quantify their oscillations. From observations, increasingly detailed knowledge of thermodynamical 
properties in filaments is available. This relies on inverting spectroscopic data, where departures from local thcrmody- 
£3 , namic equilibr ium in the radiative t ransfer through one- and two-dimensional promine nce models are n owadays routinely 
incorporated (|Labrosse et alj [20101) . At the same time, high resolution observations dLin et al.l 20051). compu tations of 



X 



three-dimensional (3D) magnetohydrostatic equilibrium configurati ons without gravity ( Aulanier et al.l l 2002D , or mod 



ern magnetic field extrapolations fr om photospheri c magnetograms dDeRosa et al.l[2009l ). together with the application 
of spectro-polarimetric techniques (|Lin et al.l 119981 : iTruiillo Bueno et al.l 12002^ . have brought insight into the complex 
three-dimensional topology of the surrounding magnetic field. Simultaneously, increased spatio-temporal resolution has 
quantified a host of oscillation modes observed in prominences, with velocity amplitudes of a few km s . Both short- and 
long-period signals, from minutes to hours, can be detected, and regions that oscillate cohere ntly with similar periods 
allow us to quantify damping rates for seemingly global eigenoscillations (|Terradas et al.|[2002T) . 

Theoretical studies in prominence seismology have meanwhile identified the most likely candidates explaining these 
motions, in terms of fundamental sausage, kink, and other mode types present in the MHD spectrum of magnetized slabs 
or cylinders. To explain the observed damping of the eigenmodes, one can on the one hand study wave damping as a 
result of non-adiabatic effects, such as anisotropic, thermal conduction, and radiative losses. For cylindrically structured 
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coronal loops and prominences, with prevailing helical fields, early work bv lKeppens et all (|1993|) demonstrate d how these 
non-a diabatic effects can give rise to damped and/or overstable global, discrete non-adiabatic Alfven waves. ISoler et all 
(|2008t) revisited these effects for cigenmodes of more simplified, uniform cylinder models in homogeneous magnetic fields. 
By using uniform internal/external conditions, the possibility of wave resonant behavior was excluded artificially. This 
resonant behavior, where global waves have eigenfrequencies in the range of the MHD continuous spectrum, provides an 
alternative means of explaining t he observed damping rates. The theory for resonant MHD waves is revi ewed in many 
pape r s, e.g. Goossens et al.l (1201 if ), and references therein. In other reviews of prominence oscillations fe.g. lMackav et all 
l2010HArregui fc Ballesterll201ir i, many of the accumulated theoretical findings are discussed, and a striking observation 
is that all these efforts dramatically simplify the underlying magnetic topology, in order to allow analytic progress. Non- 
uniformity, and t he associated possibility of resonant absorption, is at most introduced by varying the density profiles (e.g. 
ISoler et aLll2010l ). sticking to uniform field topologies. At the same time, efforts to ac count for the mult i-stranded nature 

of actual solar filaments have also been initiated, in similarly simple magnetic fields dLuna et al. 20101). 

To pave the way for more detailed prominence seismology, true MHD spectroscopy ( Goedbloed et al.lll993l ). of the 
kind pioneered in fusion-related plasma physics, is called for. This relies on accurate computations of the underlying 
equilibria, with full account of the complex magnetic topology, and subsequent computation of its eigenspectrum. In 
our accompanying paper I, we started this effort by demonstrating our ability to compute accurate translationally 
symmetric flux rope equilibria containing prominence condensations, governed by the magnetohydrostationary equations 
balancing Lorentz forces, pressure gradients, and gravity. We here continue this effort, and use these equilibria, where 
we incorporate helical fields, axial shear flows, and varying plasma beta conditions consistent with observed filament 
properties. We here concentrate on quantifying all eigenmodes defining the continuous parts of the MHD eigenspectrum. 
This is a worthwhile effort in its own right, partly because the global, discrete modes are known to be organized about 
this continuous spectrum (see e.g. modern textbooks such as iGoedbloed fc Poedtsl 120041 : iGoedbloed et al.ll2010D . and 
because the possible damping of these global modes due to resonant absorption relies on knowledge of the frequency 
ranges that they occupy. 

Charting out the continuous MHD spectrum for realistic flux rope structures with prominences has never been 
attempted before, in sharp contrast to the vast body of knowledge on MHD eigenmodes in axisymmctric, tokamak 
configurations. Our effort considers flux ropes idealized to be invariant along the filament axis, with an outer cylindrical 
cross-section, for simplicity. Since the gravitational stratification will act to deform a pure cylindrical flux model to one 
with still nested, but downwards displaced flux surfaces, the resulting poloidal variation in the equilibrium quantities will 
induce coupling between the different branches of the continuous MHD spectrum. These branches are labeled by their 
axial, k, and azimuthal, to, mode numbers. When charting the mode frequencies from the filament axis to the flux rope 
edge, the spectrum can help us to determine the avoided crossings, which are akin to the poloidal m ode couplings discussed 
for equilibria without gravity, such as the nested elliptic cylinders studied bv lChan.ee et al.l (119771). or configu rations with 
non-circular cross-sections and/or density variations along field lines analyzed bv lPoedts fc Goossens! (jl99lh . We extend 
these latter studies to translationally invariant, gravitating equilibria with possibly shear flow v z (i)j), restricted to be a flux 
function. The actual equilibria are presented and discussed in paper I. The analysis of their continuous spectr a presented 
here, benefits from the detailed knowledge of the spectrum of the more complex, axisymmetric t opologies (iGoedbloedl 
19751 ). The case for self-g ravitating, static axisymmetric magnetized tori was treated analytically in lPoedts et al.l (|1985[ ). 
Goedbloed et al.l (|2004bl ) analyzed of the continuous MHD spectrum for magnetized accretion tori, where external gravity, 



nested toroidal flux surfaces, and both poloidal and toroidal (Keplcrian-like) rotation is in cluded, and identified a p otential 
for unstable localized modes called trans-slow Alfven continuum modes. The approach of lGoedbloed et al~l (|2004bt ). which 
combined analytic findings for the equilibria, the poloidal mode couplings that were confronted with numerically computed 
equilibria, and eigenmodes is analogous to our efforts in this series. 

Our paper is then organized as follows. In Sect. [21 we present the equations governing all eigenmodes for these 
gravi tating, 2.5D MHD equilibria. Since we allow for axial flow, the Frieman-Rotenberg formalism (Fricman fc Rotenberd 
ll960T) is revisited, and using the straight field line coordinates (see paper I) describing our equilibria, we present them in a 
form where a field line projection can be made. This form can directly be compared with all other cases studied previously, 
such as the cylindrical limit, or the analogy with tokamak configurations that contain toroidal flow. In Section[3J we then 
derive the equations for the continuous spectrum, and find that prominence equilibria where density is a fl ux function can 
also be affected by convective continuum instabilities dBlokland et al.ll2007at Ivan der Hoist et ai]|2000bl lal). The coupled 
system of differential equations can be converted to algebraic relations between poloidal Fourier harmonics, which clearly 
identifies how Alfven and slow modes can couple. This coupling is then studied analytically in Section |4j where the small 
gravity parameter expansion previously used (paper I) to analyze equilibrium effects due to gravity, re-enters to analyze 
the leading order coupling effects. For low frequencies, we then demonstrate that a gap can arise in the Alfven continuum 
branches, from a four-mode coupling scheme. In Section [5j we briefly summarize the numerical approach used in the 
remaining sections. There, the equilibria presented in paper I are all diagnosed spectroscopically, by charting out their 
continuous spectrum. The analytic findings are used to check and analyze in detail the derived eigenfrequencies. This is 
done for cool prominences embedded in hot cavities in section [6j and for complex, double layered prominences without 
and with axial flows in section [7] Conclusions and an outlook to future work is given in section [5] 



2. Spectral formulation 



In what follows, we analyze the equations governing linear, ideal MHD perturbations about a stationary equilibrium. The 
equations for the equilibrium have been presented and solved numerically in paper I. The (time-independent) equilibrium 
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quantities are denoted by the magnetic field vector B, the pressure p, the density p, and the external gravitational potential 
<I>, which for the prominence application at hand are all functions of the poloidal coordinates The background 

equilibrium may contain a purely longitudinal shear flow v = v z (ip)e z , where the axial flow is constant on each flux 
surface. Each flux label ip identifies a single flux surface, which form nested surfaces defined by the magnetic field 
topology within the prescribed outer flux rope shape. We first recapitulate the governing equations, and then present a 
representation exploiting a projection involving straight field line coordinates. 

2.1. Frieman-Rotenberg formalism 

The formalism introduced bv lFrieman fc Rotenberel ([i960) is used here to investigate the stability properties of stationary 
MHD equilibria. From the full set of linearized MHD equations, these authors derived one second-order differential 
equation for the Lagrangian displacement field vector £ 

<9£ <9 2 £ 

G(e)-2pvvJ-p^ = 0, (1) 

where the generalized force operator G(£) is 

G(£) = F(£) + V$V • 0£) + V ■ [p£v • Vv - pvv • V£] . (2) 

Here, 

F(£) = -Vn + B-VQ + Q-VB (3) 

is the force operator for static equilibria without gravity, which was derived bv lBernstein et all (|l958l ). These expressions 
contain the Eulcrian perturbation of the total pressure 

n = - 7P V • £ - £ ■ Vp + B • Q, (4) 
where 7 is the ratio of specific heats, and the Eulerian perturbation of the magnetic field 

Q = V x (£ x B) . (5) 

The time-dependence of the displacement field is assumed to be exponential with normal mode frequencies u>, £ = 
£exp(— kot). Using this assumption, the Frieman-Rotenberg equation can be written as 

G(£)+2ipwv- V£ + pcj 2 £ = 0. (6) 

Wc base our derivations below on the axial symmetry of the solar prominence equilibrium by writing the axial dependence 
of the displacement field as £ ~ exp(ifcz), where k is the axial mode number. We then neglect the hat on £ and suppress 
the implied exp(— icut + ikz) dependence for all linear quantities. 



2.2. Field-line projection and representation 

As in all stability studies in fusion research, we need to exploit a projection based on the nested topology of the magnetic 
surfaces, on which the helical field lines reside. This projection ensures that on each flux surface, one can make a 
distinction between two wave directions: one parallel and the other perpendicular to the magnetic field. In the short 
wavelength limit, these directions correspond to the polarizations of the slow and Alfven waves, respectively. To perform 
this pr ojection, we adopt straight fi eld line, non-orthogonal coordinates (i/j, d, z) that are introduced in our accompanying 
paper (jBlokland fc Keppend[20Tl . which are known once the equilibrium is fully determined. The triad of unit vectors 
based on the projection is 

Wib B 
n "^' ^ bxn < b ^s< (7) 

where the gradient of the magnetic flux ip is perpendicular to the flux surfaces. Using these unit vectors, the components 
of the displacement field £ can be replaced by the three related variables 

X = B*e-n, r^i^-£-7r, Z = i-U-b, (8) 



and the projected gradient operators become 



2?=-^n-V =oV-— (9) 



Q = —IBftBTt ■ V = —j-d# — kB$, (10) 
^=-iB-V =^d<5 + ^- (11) 
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The second equality in the expressions for JF and Q only holds when these operators act on a component of the displace- 
ment field These expressions contain several quantities already introduced in paper I, such as the metric elements gij 
and the Jacobian J for the non-orthogonal straight field line coordinates, and quantities shown to be flux functions to 
realize equilibrium, such as the poloidal stream function / — B z and the safety factor q(ijj) = 31. 
The Frieman-Rotenberg equation in Eq. (TTJ) can be written as 



(A + G)x — pu; 2 Bx = 0, with x 



(12) 



by exploiting the straight field line coordinates and applying the projection defined by Eq. ([5]). where u(tp) = w — kv z (ip) 
is the flux-surface localized Dopplcr-shiftcd cigcnfrcquency. Furthermore, the matrices A and B arc 3x3 matrix operators 
that are also present in the case of a static equilibrium without gravity. The matrix operator G contains the elements 
associated with the external gravitational potential. The non- vanishing matrix elements of A and B are 



An 

Aia 
Ai 3 
A 2 i 

A 22 



-V 



IP- 



S' 



3 



V^3- 



2[V(B$k v )\ 



A 23 

A31 = T 



B 2 



G J VKJ 



2kB^K 
B 2 



p> 



(13) 



B 2- 



j 



A32 = FipG 



B 2 ' 



FjpJ 7 , 



and 



Bll = M' 



R - B # 
V22 = -Jp, 



B33 = B 



We used the poloidal curvature of the magnetic surfaces k p from paper I, along with the operators 



922/ 



OA 9 — 

922 



(14) 



(15) 



These are related to the normal gradient operator T>, originally introduced by iGoedbloedl (|1997f ) for the spectral analysis 
of static tokamak equilibria. The square brackets in the expressions in Eq. (|13[) for the matrix elements of A and the 
gradient operators defined by (fl"5"j) indicate that the differential operator only acts on the term inside the bracket. This 
notation is also used in the expressions below. As expected, the matrices A and B ca n also be de r ived f r om th e expressions 
for wave modes about static axisymmetric MHD equilibria, found in the papers by lGoedblo"edl (|l975l . Il997t) , by formally 
setting the toroidal distance R to unity. The matrix G enters if there is an external gravitational potential present. The 
expression for this matrix is 



,PV [ -A 



(win + xg) 



B 2 



{YD n + A J") p 



\ 



where 



p (ipjVU + 7 



/' 

A 



ell 

s 4 



■jdvp, 



P\ ~5t 



B 2 



ikp 























-jd^p 


1 



(16) 



1 



J 
V<$> 



d«<f> 



pJ 



-] 
P 



-B^ • Vp, 



(17) 
(18) 



We note that the term p represents the pressure variation on a flux surface. This ma trix G can also be derived from the 
expression for the matrix operator R found in the paper bv iBlokland et a l. (2007a), where axisymmetric, gravitating, 
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magnetized accretion tori with purely toroidal flow where spectroscopically analyzed. The equivalence involves (1) setting 
their radial distance R = 1, (2) producing the toroidal curvature K t = for a translationally symmetric solar prominence, 
and (3) replacing the toroidal mode number n by the axial wave numbe r k. In that analysis, as well as in the related one 
for toroidally rotating tokamak plasmas bv Ivan der Hoist et al.l (|2000aF ). another matrix C related to the Coriolis effect 
appeared, which naturally vanishes for our translationally invariant configurations. 

Ano ther interesting, related limit is the cylindrical limit known from standard textbook treat- 
ments dGoedbloed fc Poedtsl l200l . If we ignore gravity altogether, then matrix G will be zero and the operators 
T and Q will just be algebraic ex pressions. In this cas e, the spectral equation given in Eq. ([T2"j) reduces exactly to the 
well-known equation presented bv lHain fc Liistl (|1958l ). They rewrote this matrix form to one second-order differential 
equati on for the radial displacement field, assuming a ratio of specific heat s 7 = 1. This w as ge neralized by Goedbloedl 
(|l97lh for arbitrary 7, and rewritten, for cases involving equilibrium flow in lHameiril (jl98lh and lBondeson et alJ ( 1987 ). 
A cylindrical limit in volving flow and external gravity, suitable for analyzing accretion disk stability, was presented by 
iKeppens et al.l (|2002f ). and can also be easily shown to be directly related to the above expressions. 



3. Continuous MHD spectrum 

In the previous section, we derived the spectral equation in Eq. ([12]) governing all MHD waves and instabilities in trans- 
lational symmetric solar prominences. This equation is the starting point to all following MHD spectral computations. 
In particular, we can derive the equations for the continuous MHD spectrum by first considering modes localized on a 
particular flux surface ip = ipo- We can then obtain a non-singular eigenvalue problem in the local Alfven and slow-like 
cigenfunctions 77(1?), CW- This form allows us to derive a stability criterion for continuum modes, where it is found that 
for an axially symmetric prominence, gravity can cause the continuous spectrum to become unstable. By expanding the 
periodic "^-dependence of the eigenfunctions, we get an equivalent algebraic set of coupled equations. 



3.1. General formalism 



To derive the equation for the continuous spectrum, we assume that the normal derivative (d/dip) of the eigenfunctions 
is large compared to the cigenfunctions themselves and the local Doppler-shiftcd cigenfrcquency uj(ip) is finite. Using 
these assumptions, the first row of the spectral equation in Eq. f| 1 2|) can be solved approximately 



J VijX 



jp + B 2 \ B 2 



TZ\- 



1P- 



J_gy + W (± Y 

B 2 1P + B 2 \B 2 



(19) 



We note that X is small compared to either Y or Z . This implies that the continuum modes arc mainly tangential to a 
particular flux surface. Furthermore, we note that dX/dtp, Y, and Z are of the same order. 

Inserting the approximated solution given by Eq. (|19|) into the second and third row of the spectral equation in 
Eq. (|12p results in an eigenvalue problem that is independent of the normal derivative. Hence, the reduced problem 
becomes non-singular. Exploiting this property, we can write the projected displacement field components Y and Z as 



Y 
Z 



Sty 



(20) 



where S(ip — ipo) is the Dirac delta function. Here, Y and Z have been divided into an improper (^-dependence) and 
proper part (^-dependence). This kind of splitting was introduced bv lGoedbloedl (|1975l ) for static axisymmetric equilibria 
without gravity. The reduced, non-singular eigenvalue problem becomes 



(a + g - pw 2 b) 



= 0, 



(21) 



where 



B 2 JP + B 2 



I I piiB , , pi 2 „ 
- 4 ^7^ ( ^ ) + ^<P°' 



7P + B 2 
7P + B 2 



PI' 



\ B 2 7p + B 2 



B 2 + pIN 2 



m,pol 



-B 2 



1 „ PP 



/ 



7P + B 2 B 2 



JB 2 jp+B 2 



B 2 



m,pol 
2 Ar2 



pB'N^ po] 







(3. 

B 2 
1 B 2 



(22) 



(23) 



(24) 
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In the expressions above, Kg is the geodesic curvature, 



I 



JBoB 2 



(25) 



where n = b • Vb is the field line curvature. Furthermore, pol is the magnetically modified Brunt- Vaisala frequency 
projected onto a flux surface, 

pp 



N 2 , 

m,pol 



JL 

B 2 



Jp 



IP 









[ pB _ 




_ P B ' { 



VB 2 
Vp 



P 



B 2 



Vp 



(26) 



This magnet i cally mod ified Brunt- Vaisala frequency has also been found in various related cases. For example, 
iPoedts et al. (119851) and Belien et ah! (|l997t) found this frequency for static, axisymmetric gravitationall y stratified plas- 
mas. Ivan der Hoist et al.l (|2000bllaF discussed its role for tokamak plasmas with purely toroidal flow. iBlokland et al.l 
(|2007a[) generalized the expression by including an arbitrary external gravitational potential in its application to rotat- 
ing, magnetized accretion disks. When we examine the matrix g closely, we see that the poloidal variation in the pressure 
is present in all its matrix elements. This has the consequence that the diagonal dominance of matrix a will be destroyed, 
which causes possible mode coupling between Alfven and slow modes. 

The non-singular eigenvalue problem presented in Eq. (|21|) is solved for each flux surface, using poloidally periodic 
boundary conditions for the eigcnfunctions r\ and Q. Typically, one solves the problem for a given axial mode number k, 
and we then get a set of discrete eigenvalues for each flux surface. When we quantify these eigenfrequencies for all flux 
surfaces, they together map out the continuous spectra. For completeness, we point out that if one uses a formulation 
based on the primitive variables p, v, p, and B, an additional entropy continuum kv z is found. iGoedbloed et al.l (|2004aj) 
showed that this Eulerian entropy continuum, in ideal MHD, does not couple to any of the other continua, nor to the 
discrete stable and/or unstable modes. 



3.2. Spectral properties and stability criterion 



We now derive a general stability criterion for the continua. For this derivation, a Hilbert space with appropriate inner 
product needs to be constructed. The parallel gradient operator J 7 is a Hermitian operator under the inner product 



Jv* (Tw) d& = j> J {Tv) wdi9 (27) 
for poloidally periodic functions v and w. Inspired by this property, a Hilbert space can be defined with the inner product 

v* • wJdtf. (28) 



In this Hilbert space, the matrix operators a, b, and g are Hermitian operators. Setting v and w equal to the eigen- 
functions (77, C) T , one can derive a simple quadratic polynomial for the eigenfrequency ui from the spectral equation in 
Eq. ([2"TT) for the continua, namely 



aui — c = 0, 



where the coefficients are 



f|M 2 + 5 2 ICI 2 



(29) 

(30) 
(31) 



-pB 2 ^ 



2 \j2 

BV.pol 



B 2 



v + C 



Here, N^y pol is the Brunt- Vaisala frequency projected onto a flux surface 



I' 



BV,pol 

Btf • Vp 



la PP 

— dap 

Jp -ip 



pB 
B_a_Vp 

pB 



B,9 
p~B 

jBS 



Vp - — Vp 

7P 



(32) 
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The coefficient a is real and always greater than zero. As for a, the coefficient c is real but can be zero or even negative. 
The solutions of this second order polynomial in Eq. (|29[) arc simply uj = iy^c/a. If c is negative, a damped stable wave 

and an overstable mode are obtained. In this case, its absolute value is found from \uj\ = ^/—c/a, while Rc(w) = kv z , 
so the oscillation frequency is set by the local Doppler shift. By looking at the expression for c, we realize that this 
coefficient is equal to or greater than zero if 

^Iv.poi > (33) 

is satisfied everywhere in the plasma. Therefore, equilibria for which the entropy or the temperature (assuming 7 > 1) is 
a flux function, will always be stable according to this criterion. However, when density is a flux function, the equilibrium 
variation may violate this criterion. This implies that such equilibria have unstable MHD c ontinua, a situation identifi ed 
in earlier, related wo rk, as being caused by the convective continuum instability (CCI) ([van der Hoist etHl l2000blial 
iBlokland et al.ll2007al) . This CCI instability criterion in Eq. (|33p is analogous to the Schwarzschild criterion for convective 
instability, with one important difference. The Schwarzschild criterion deals with nor mal derivatives, while thi s criterion 
deals with tangentia l ones . This difference has also been found by vari o us aut hors (jHellsten fc Spies! (|1979( ): lHameiril 
(jl983h ; |Poedts et al.l (|1985f ); Ivan der Hoist et all (|2000al ) ; iBbkland et all (|2007aD ). 

The relevance of this CCI instability to solar prominences may well be significant. We note that it is common practice 
to study the (slow) evolution of prominences carrying flux ropes by taking an equilibrium configuration, and subjecting 
it to e.g. photosph eric-boundary-driven changes, and relax it thereby repeatedly to successive equilibrium states (see e.g. 
lAmari et al.l ll999l). As a result, the internal flux rope variation proceeds through a series of equilibria. For short-lived 
prominences, equilibria with density as a flux function may occur, and in this evolutionary process, the equilibrium can 
transit from the stable to the unstable regime as governed by the CCI criterion. As we later demonstrate, the CCI growth 
rate also varies with increasing gravity parameter g. As explained in Paper I, this p-increase could relate to an increase 
in size, internal density, or the weaker overall magnetic field strength of the prominence carrying flux rope. 



3.3. Expansion in poloidal harmonics 

We now continue to transform the two coupled differential equations for the MHD continua in Eq. (f2~Tj) to an algebraic 
set, by expanding the eigenfunctions rj and £ in poloidal Fourier harmonics. This expansion has the advantage that we 
can perform a detailed investigation of how different poloidal harmonics may couple to each other. This coupling could 
cither be between modes of the same kind, for example Alfven-Alfvcn coupling, or between modes of different kind such 
as Alfven-slow coupling. In the next section, we show that in an actual solar prominence, more complex mode coupling 
is present. 

Expanding the eigenfunctions 77 and C for the MHD continua in a finite Fourier series, 



E 



(34) 



and applying the operator 



Je 



-imi? 



do 



to the spectral equation in Eq. (|21[) . one obtains the algebraic coupled equations 

= 0. 



L, 



(35) 



(36) 



The operator in Eq. (1351) is the result of exploiting the inner product in Eq. (|28[) of the previously defined Hilbert space. 
The poloidal mode numbers m and ml have the same range and the 2x2 submatrices are given by 



E 



Li, 



(37) 



The explicit expressions for the integrand functions K{{f), £($), M($), and JV(i?) are 
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These expressions are derived by exploiting the Hermitian property in Eq. (|27[) of the parallel gradient operator J 7 . 

Once again, we consider the cylindrical limit in combination with zero gravity, which is appropriate for many earlier 
seismological studies of prominences. In this limit, all equilibrium quantities are independent of the poloidal angle 
therefore the geodesic curvature will be zero. This means that the Alfven and slow continuum decouple fully from each 
other. Furthermore, only one single harmonic is required for the expansion of the Fourier series of the eigenfunctions r\ 
and £. Solving the spectral equation for this case, one finds the well-known local Dopplcr-shifted Alfven continua 

and the local Doppler-shifted slow continua 



^ s ,m = kv z ±u> SiTn , 

Looking at the expressions in Eqs. ([3"51) - (|4"Tj) for the submatrix elements, we can make some general statements about 
the MHD continua of solar prominences. In the absence of gravity (ji = pol = 0), the Alfven and slow continuum can 
be coupled by the geodesic curvature present in the off-diagonal terms M($) and iV($). With gravity, these off-diagonal 
terms can become more pronounced. The second statement is that an Alfven continuum at a resonant surface, associated 
with the element K($), may shift away from the local Dopplcr-shifted frequency kv z because of the (1) geodesic curvature, 
involving its square K g , (2) compressibility and finite pressure (72?), and (3) the pressure variation on a flux surface. All 
these effects suggest that an Alfven gap can appear around the local Doppler-shifted frequency and that Alfven-slow 
coupling should be taken into account for these gaps. This kind of coupling is investigated in the next section. 

4. Small gravity expansion 

In previous sections, we have derived equations for the continuous MHD spectrum of a translationally symmetric equi- 
librium. These derivations were completely general, and allowed us to identify an instability criterion on the one hand, 
as well as to identify the terms responsible for Alfven-slow type coupling. In this section, we quantify this coupling effect, 
using the small gravity parameter expansion for the equilibrium, similar to those discussed in paper I. In the remainder 
of this paper, we assume a gravitational potential in which the prominence is embedded given by 

§(«,») = {x-x )g, (44) 

where xq is the location of the center of the last closed flux surface of the prominence and the role of gravity is represented 
by the constant g. This gravity parameter g denotes a dimensionless quantification of the relative importance of gravity 
on the equilibrium structure, and we explained in paper I that this parameter can meaningfully vary by orders of 
magnitude for solar prominence configurations. For small g and assuming a circular cross-section of the prominence 
plasma, we quantified in paper I the change induced in the equilibrium (the downward shift of the flux surfaces), and 
the corresponding straight field line coordinates. This kind of expansion is similar to the small inverse aspect ratio 
expansion for tokamak plasmas. Using this analytical result for the equilibrium, we can now demonstrate that gaps (or 
even instabilities) may appear in the continuous MHD spectrum because of mode coupling, as a direct consequence of 
the force of gravity. 

4.1. Continuous MHD spectrum for equilibria weakly affected by gravity 

We apply the small gravity expansion to the expressions found in Eqs. (j3"5)) . (|39|) . (|4T)|) . and (|4I[) for the matrix ele- 
ments K($), M($), iV($), and L{&), respectively. For this expansion, we rely on the metric elements that are valid 
for straight field line coord inates exploiting the Shafranov shift A(r), given by Eq. (26) in our accompanying paper 
(|Blokland fc Keppensl[2rjTlT) . In addition, one can derive the quantities 
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where the total magnetic field Bq = ip' 2 + I 2 . The parameter t g = 1 has been introduced to track the linear and quadratic 
dependence on the geodesic curvature. Furthermore, the parameter t e allows us to distinguish the three options for the 
flux function dependence in the equilibrium, namely 



r l for TO), 
for p(tp), 
.1/7 for Sty). 



(51) 



These flux functions lead to different varia nts of the extended Grad-Sh afranov equation in Eq. (10) for the equilibrium, 
as presented in our accompanying paper (iBlokland fc Keppensl l20l"Tt ) . We again note that the magnetically modified 
Brunt- Vaisala frequency pol from Eq. (|50|) can only be negative if the density is a flux function. This once more 
confirms that for this particular choice of the flux function, the MHD continua may become unstable as already found 
to be true in general above. 

Applying the small gravity expansion, the eigenvalue equation in Eq. ((36|) can be written as 



0. 



(52) 
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We note that the first and second rows of the eigenvalue equation in Eq. ([52]) have been multiplied by B 2 /rI and q^ 1 , 
respectively, to get rid of the dimensional factors in front of the square brackets in the definitions given in Eqs. (|38p - (|4"Tj) . 
The approximated expressions of the integrands of the matrix elements used in the spectral equation in Eq. (f5"2"j) are now 
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Once more, the magnetically modified Brunt- Vaisala frequency pN^ , (jBU)) is only negative if the density is a flux 
function (t e = 0). 
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4.2. Low frequency Am = gap 



The approximate relations just given can now be analyzed in terms of mode coupling possibilities. We focus on low 
frequencies, ui = O(g), and show that a Am = gap can be created in the neighborhood of a rational surface m + kq = 0, 
where q quantifies the flux-surface specific local twist in the field lines on this surface. For such a surface, we determine 
the leading order coupling on an Alfven continuum mode rj m with the same poloidal mode number m. Close to the 
rational surface, we can safety neglect the third term of the expression in Eq. ((55|) for K and Eq. (|58|) for L. We are 
left with the eigenvalue equation for the MHD continua in Eq. (f5^j) . which reduces in essence to a four-mode coupling 
between the central Alfven mode, and the three neighboring slow modes given by 
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In this derivation, we made use of the finding that the integral of sin($) and sin 2 ($) over d gives \\(5 mim i+i 
and -j(S m , m '-2 - 2S„ 
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+ 5 m ,m'+2), respectively. We note that for this four-mode coupling scheme, we neglected any 
coupling to the mil Alfven modes, an assumption that is valid as long as the associated frequencies are far away in the 
spectrum. The same argument also holds for the coupling to the m ± 2 Alfven and slow modes. By careful examination 
of this 4x4 matrix, one notices that the coupling between C m ±i W1 th either rj m or £ TO is on the order of O(g), while all 
the other types of coupling are on the order of 0(g 2 ). Because the main modes involved have the same poloidal mode 
number m, this coupling is known as a Am = gap possibility. 

The matrix of the reduced eigenvalue problem presented in Eq. (|61[) is symmetric and therefore its eigenvalues uj 2 
will be real. This is consistent with the properties of the quadratic polynomial (|29[) for the eigenfrequency ui. If gravity 
drives a certain MHD continuum unstable, it means that the corresponding eigenfrequency uj 2 will be smaller than zero, 
resulting in a complex eigenvalue u;. Once the eigenvalue uj is known, it is straightforward to determine the dominant 
character of the MHD mode by determining the dominant component of the eigenvector of the eigenvalue equation in 
Eq. (|6ip . This analytic result will be used in what follows to verify the numerically obtained continua. 



5. The spectral code PHOENIX 

While the equilibria presented in paper I were computed with the FINESSE code (jBelien et al. 20021). the stabilit y 
analysis of the prominence models is subsequently done using the spectral code PHOENIX ( Blokland et al.l l2007bh . 
Although this code was originally developed for analyzing axisymmctric plasma tori such as tokamaks and accretion 
disks, the design is flexible enough to analyze axially symmetric, gravitating plasma tubes. This possibility is exploited 
in the present PHOENIX version 1.2. For a given two-dimensional equilibrium computed by FINESSE and converted to 
straight field line coordinates, we then need to select an axial wave number k and a range of poloidal wave numbers m 
for which we determine the complete MHD spectrum numerically. A range of poloidal wave numbers is required to allow 
for the anticipated mode coupling caused by the combinations of non-circular cross-sections of the flux surfaces and/or 
the presence of a gravitational potential. Instead of solving the spectral equation in Eq. (fT3i)) , PHOENIX solves the 
linearized version of the full set of non-linear MHD equations. This allows us to extend the cigenspectrum computations 
for future incorporation of non-adiabatic effects, while it in addition will contain the trivial solution for the entropy 
modes known as the Eulerian entropy continuum. The discretization of these linearized MHD equations exploits a 
finite element method in the radial direction, while in the poloidal direction a spectral method is applied. The finite 
elements used are a combination o f quadratic and cubic Hermite elements. This combination prevents the creation of 
spurious eigenvalues (|Rappazl 11977). The resulting generalize d non-Hcrmitian eigenvalue problem is then solved using 
the Jacobi-Davidson method ( Sleiipen fc van der Vorstl Il996l) . For such an eigenspectrum computation, one needs to 
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supply boundary conditions at both the magnetic axis (set by symmetry) and the outer flux surface. There, we assume 
a perfectly conducting wall for simplicity. Strictly speaking, this condition is not applicable to true flux ropes carrying 
embedded solar prominences. However, since we focus on the continuous spectrum alone, the boundary conditions in the 
radial direction have hardly a ny influence, since t h ose w aves or instabilities are by definition localized on flux surfaces. 

Exploiting this realization, iPoedts fc Schwarzl (119931 ) introduced a more efficient method than the Jacobi-Davidson 
method to compute the continuous MHD spectrum. This method replaces the radial dependence of the cubic Hermite 
and the quadratic elements with log(e) and 1/e, respec tively, on each ind ivid ual flux sur face. In this way, the singular 
behavior of the continuous modes can be approximated. iGoedbloedl (|l975h and lPaol (ll975h describe in d e tail the reasons 
for this type of singular behavior in the continua for static plasma tori without gravitv TGoedbloed et al.l (|2004b[ ) pointed 
out that similarly singular behavior for the MHD continua prevails in plasma tori including toroidal plus poloidal flows 
and gravity. The same behavior is thus also expected for our translationally symmetric case. After replacing the elements, 
the eigenvalue problem becomes a small algebraic problem on each flux surface, with a size determined by the number 
of poloidal mode numbers taken along. For each flux surface, this can be solved using a direct QR method. All spectra 
below are computed in this way, for an axial wavenumber k = — 1 and using seven poloidal mode numbers centered 
around m = 1. A larger number of poloidal Fourier harmonics can be used, which results in a more accurate solution 
but also in a longer computation time. The typical radial resolution used is either 101 or 1001 points, depending on how 
much detail is required in identifying each branch of the MHD continua. 



6. Cool prominence surrounded by hot medium 

The theory presented will be applied to the two equilibrium classes presented in paper I, cool prominences embedded in 
hot medium and double layered prominence s. We now discuss the first class. The equilibria and their associated profiles 
are discussed in our accompanying paper (jBlokland &: Keppensl [201 it ). As described in our accompanying paper, we 
analyze the continuous MHD spectrum for increasing values of the gravity parameter g. 



6.1. Gap creation when increasing g 

Fig. [T] shows part of the continuous MHD spectrum as a function of the radial flux coordinate s = \p§ computed by 
PHOENIX for a prominence equilibrium with gravity g = 0.001 and assuming that the temperature is a flux function. The 
frequencies are normalized to the Alfven time ta = L^JJIq p mag / B mag , where [1q, L, /9 mag , and -B mag are the permeability, 
the radius of the prominence, the density on the magnetic axis, and the magnetic field on the magnetic axis, respectively. 
Near the magnetic axis where s = 0, each branch can be labeled immediately using the analytical expressions in Eqs. (|42l) 
and (|43|) for the MHD Alfven and slow continua. At the location s w 0.9, the safety factor q of the equilibrium becomes 
the value of one (see paper I). At this location, the m = 1 Alfven and the m = 1 slow continuum all become zero. At the 
low gravity value for Fig. [I] no true gap can be detected, and for progressively stronger gravity, we expect mode coupling 
to occur, ultimately inducing a possible gap or even instabilities in the continua. At the same location, the to = slow 
and to = 2 slow continuum have the same frequency, namely approximately 0.797. 

We next take a closer look at the continuous MHD spectrum about the q = 1 surface, to compare the PHOENIX 
results with the ones predicted by the four- mode coupling scheme in Eq. (|6ip . This comparison is presented in Fig. [2] 
Results from PHOENIX are represented by crosses, while the colored lines result from solving the eigenvalue equation in 
Eq. (|6ip . The colors indicate the dominant characters, which are easily determined from the associated eigenvector. The 
colors red, green, blue, and yellow correspond to dominant rj 1 (Alfven), (slow), £ (slow), and £ 2 (slow) characters, 
respectively. The plot shows that there is excellent agreement between both methods, as expected for small gravity. 
Furthermore, the m = 1 Alfven continuum (red) intersects the to = (blue) and to = 2 slow continuum (yellow). This 
intersection may lead to mode coupling, resulting in a Am = 1 gap. The analysis of this kind of gap is left to future work. 
The crosses at zero frequency correspond to the Eulerian entropy continuum, as also found numerically by PHOENIX. 

We now repeat the same analysis, for a case with intermediate gravity g = 0.500 and still assuming that the temper- 
ature is a flux function. The results of PHOENIX results are in excellent agreement with the solutions of the four-mode 
coupling scheme in Eq, (foTj) . as shown in Fig. [3] For this larger g value, the results now clearly illustrate that a gap 
appears in the m = 1 Alfven continuum, while no gap appears in the to = 1 slow continuum. We note that the radial 
flux coordinate s of the q = 1 surface is shifted inwards relative to the previously discussed equilibrium. However, the 
radial location r of this surface remains almost the same (see paper I). 

If we continue to increase the gravity parameter g, for a cool prominence equilibrium where temperature is a flux 
function, we find progressive deviations between numerical and (approximate) analytical predictions. The plot presented 
in Fig. U shows the results from PHOENIX (crosses) and the four-mode coupling scheme (colored lines). Near q = 1, 
the results of the four-mode coupling scheme corresponding to the Alfven continuum (red) show good agreement with 
PHOENIX, while for a slow continuum (green) there is a discrepancy. The four-mode coupling scheme indicates that 
a gap in the slow continuum may result, while PHOENIX does not find such a gap. The gap in the Alfven continuum 
clearly becomes larger if gravity is increased. The radial location r of the q = 1 surface remains almost the same as in 
the two previously discussed cases, while its radial flux coordinate s is shifted outward. 

In the previous three equilibria, which only differ in terms of the increased role of gravity, we observed that the 
gap appearing in the Alfven continuum becomes larger if the gravity is increased. This dependence was investigated 
by determining the Alfven frequency at the q = 1 surface, for various values of the gravity parameter. The result of 
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Fig. 1. Part of the continuous MHD spectrum for a cool solar prominence with gravity parameter g = 0.001 computed 
by PHOENIX. Shown are eigenfrequencies normalized to the Alfven time ta as a function of the radial flux coordinate 
s = all corresponding to flux surface localized eigenoscillations. 




0.80 0.85 0.90 0.95 1.00 

s 

Fig. 2. The real part of the sub-spectrum of the MHD continua of a cool solar prominence equilibrium with small gravity 
(g = 0.001) as a function of the radial flux coordinate s = for axial wavenumber k = —1 and poloidal mode numbers 
m = {—2, ...,4}. The temperature is assumed to be a flux function. The crosses are PHOENIX results, while the colored 
lines compare them with the four-mode coupling scheme. The colors red, green, blue, and yellow correspond to dominant 
771 3 Cij Co) an d C2 components of the eigenvector, respectively. 



this parametric investigation (involving the construction of numerical equilibria, subsequently diagnosed for their MHD 
continua by PHOENIX) has been plotted in Fig. [5j The plot shows excellent agreement between PHOENIX and the 
growth of the 'gap' predicted by the four-mode coupling scheme. It is seen that the Alfven gap frequency depends 
non-linear ly on the gravity parameter g. Despite this excellent agreement for the Alfven continuum gap up to g = 1, 
successive differences occur, e.g. as noted for the slow continua shown in Fig. [4] In Fig. [5j the horizontal axis does not 
use the gravity parameter g itself, but rather a scaled value. The equilibrium with gravity g = 0.001, g — 0.500, and 
g = 1.000 correspond to rp g/p = 0.0003, rp g/p = 0.154, and rp g/po = 0.300, respectively. 
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Fig. 3. The real part of the sub-spectrum of the MHD continua of a cool solar prominence equilibrium with intermediate 
gravity g = 0.500 as a function of the radial flux coordinate s = y/rp for axial wavenumber k = — 1 and poloidal mode 
numbers m = {—2, ...,4}. Temperature is assumed to be a flux function. Crosses are the PHOENIX results, while the 
colored lines give the four-mode coupling scheme. The colors red, green, blue, and yellow correspond to the dominant 
rj l , C 1; C , and C 2 components of the eigenvector. 
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Fig. 4. The real part of the sub-spectrum of the MHD continua of a cool solar prominence equilibrium with strong gravity 
(g = 1.000) as a function of the radial flux coordinate s = s/ip are shown for axial wavenumber k = — 1 and poloidal 
mode numbers m = {—2, ...,4}. Here, the temperature is assumed to be a flux function. The crosses are the PHOENIX 
results, while the colored lines are the results of the four-mode coupling scheme. The color red, green, blue, and yellow 
correspond to the dominant rj 1 , C 1; £ , and C2 components of the eigenvector. 



6.2. Unstable continua 

Finally, to illustrate the effect of the chosen flux function in the equilibrium on the MHD continua, we now discuss in 
the intermediate gravity case g — 0.500, the continuous spectra when the density or the entropy is a flux function. We 
recall that when the density is a flux function the continuous MHD spectrum may become unstable. The continuous 
MHD spectrum of the density case is presented in Fig. [5] The figure shows there is again an excellent agreement between 
the four- mode coupling scheme and PHOENIX. The plot shows two distinct features, that are not present when the 
temperature is a flux function. Firstly, no gap appears in the m = 1 Alfven continuum. Near the q = 1 surface (s ~ 0.90), 
an interesting observation can be made, where Alfven and slow mode characters mingle. A mode at for example s = 0.89 
with dominant m — 1 Alfven character slowly changes to a dominant m = 1 slow character when it is traced closer 
to the q = 1 surface. The transition points are indicated by the arrows in the plot. The second distinct feature is the 
merger of the two m = 1 slow branches to one at s w 0.852 and s « 0.94. Exactly at these locations, the m = 1 slow 
continuum turns unstable as shown in Fig. [71 This figure quantifies the imaginary parts of the continuum eigenfrequencies, 
co rresponding to (normal ized) growth rates. This instability is the convective continuum instability (CCI) as discussed 
bv lBlokland et al.l (|2007al) . There is again excellent agreement between PHOENIX and the four-mode coupling scheme. 
Furthermore, the transition points of the dominant character of the mode indicated by arrows in Fig. [7J are exactly at 
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Fig. 5. The frequency of the Alfven continuum mode rj 1 at the resonant q = 1 surface as a function of the gravity. The 
temperature is assumed to be a flux function. 
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Fig. 6. The real part of the sub-spectrum of the MHD continua of a cool solar prominence equilibrium with intermediate 
gravity [g = 0.500) as a function of the radial flux coordinate s = ^/ip are shown for axial wavenumber k = — 1 and 
poloidal mode numbers m = {—2, 4}. Density is assumed to be a flux function. The crosses arc the PHOENIX results, 
while the colored lines are the results of the four-mode coupling scheme. The color red, green, blue, and yellow correspond 
to the dominant fj x , C l7 C , and ( 2 component of the eigenvector. 



the same locations as those indicated in Fig. [SJ From this analysis, we conclude that the CCI at the q = 1 surface has a 
dominant Alfven character. 

As in the case where temperature is a flux function, we can parametrically investigate the dependence of the growth 
rate of the CCI on the gravity parameter g. This analysis again involves repeated FINESSE computations for increasing 
values of g, all diagnosed with PHOENIX, and the extraction of the growth rate for the most unstable mode in each case. 
These results are presented in Fig. [5] Once more, there is excellent agreement between the two methods. Furthermore, our 
analysis clearly identities a linear dependence between the growth rate and the gravity parameter g. Owing to our non- 
dimensionalization, this implies that there is a larger growth of the CCI mode for flux ropes with embedded prominences 
that (1) gradually grow in cross-sectional extent, and/or (2) become more heavy (density increase), and/or (3) decrease 
in magnetic field strength. 

In the context of prominence dynamics, the existence of equilibria via ble to unstable con tinua may have relevance 
to the complex dynamics observed in high resolution prominence movies (|Berger et al.ll2008l ), as well as to the sudden 
disappearance of prominence plasma, the so-called disparition-brusque phenomenon (see e.g. lTandberg-Hanssenll201ll . 
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Fig. 7. The imaginary part of the sub-spectrum of the MHD continua of a cool solar prominence equilibrium with 
intermediate gravity (g = 0.500) as a function of the radial flux coordinate s = s/ip are shown for axial wavenumbcr 
k = — 1 and poloidal mode numbers m = {—2, ...,4}. These quantify the growth rates of the CCI instabilities. Density is 
assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the resulte of the four- 
mode coupling scheme. The color red, green, blue, and yellow correspond to the dominant rj 1 , £ 1; £ , and £ 2 component 
of the eigenvector. 



0.20 - 



0.15 - 



lm(£0T A ) 



0.10 



0.05 



0.00 



- 1 1 1 — 



Growth rate of Convective Continuum Instability 



- 1 — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — 



— four mode coupling scheme 
+ PHOENIX 



_l I I 



_l I I I I I I I I I I I I I I I I I I I I I 



0.00 0.05 0.10 0.15 0.20 

r Po9 / Po 



0.25 



0.30 



Fig. 8. The growth rate of the CCI mode at the resonant q = 1 surface as a function of the gravity parameter. The 
density is assumed to be a flux function. 



and references therein). It should be emphasized that the onset of the CCI instability will be localized on those flux 
surfaces where the CCI criterion is violated. As shown in Fig. [JJ growth rates can vary from one flux surface to another 
in a given unstable case, while in addition, the maximal growth rate also varies because of gradual equilibrium changes, 
as demonstrated in Fig. [5] Whether this linear prediction for the onset of internal, continuum-driven dynamics leads to 
the rather turbulent prominence dynamics revealed by contemporary observations, or even to full disruption and the 
disappearance of prominence structures, must be answered by follow-up high resolution, nonlinear numerical simulations. 
This requires initializing multi-dimensional nonlinear simulations with the spectrally diagnosed equilibria. 

6.3. Spectra for equilibria in which entropy is a flux function 

Wc now consider the case where entropy is a flux function, for an equilibrium with intermediate gravity (g = 0.500). The 
resulting continuous MHD spectrum is presented in Fig. |H1 There is again excellent agreement between both methods. 
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It can be clearly seen that no gap in the Alfven continuum appears. The continua for the strong gravity g = 1.000 case 
have also been computed, which also show no appearance of a gap. This can be understood by examining Eq. (|29[) for 
the eigenfrcquency. When the entropy is a flux function, the Brunt- Valsala frequency projected on a flux surface -/Vbv.poI 
is zero, thereby transforming the eigenfrcquency given by Eq. (|29[) . to one similar to a case without gravity, for which no 
gap can appear. 
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s 

Fig. 9. The real part of the sub-spectrum of the MHD continua of a cool solar prominence equilibrium with intermediate 
gravity (g = 0.500) as a function of the radial flux coordinate s = yftfi are shown for axial wavenumber k = — 1 and 
poloidal mode numbers m = {—2, ...,4}. Here, entropy is assumed to be a flux function. The crosses are the PHOENIX 
results, while the colored lines a£e the results of the four-mode coupling scheme. The color red, green, blue, and yellow 
correspond to the dominant rj 1 , C 1 , ( , and £ 2 component of the eigenvector. 



7. Double layered prominences 

We discuss the continuous MHD spectrum of the second equilibrium class, the double layered pro minences. As for cool 
prom inences, the equilibria and the associated profiles are discussed in our accompanying paper (|Blokland Sz Keppensl 
l201lh . For this equilibrium class, we also investigate the influence of the gravity on the continuous MHD spectrum by 
slowly increasing the gravity. 

7.1. Continua for increasing gravity 

The first equilibrium that was analyzed with PHOENIX is the one where the temperature is assumed to be a flux function 
and the gravity is low at g = 0.001. Part of the MHD spectrum is shown in Fig. [TUJ The individual branches near the 
magnetic axis can be labeled using the analytical expressions for the Alfven continuum in Eq. (|42[) and slow continuum 
in Eq. (|43|) . The associated radial flux coordinate s of the q = 1 surface is approximately 0.927. At exactly this location, 
the m = 1 Alfven and m = 1 slow continuum are zero, as expected, and the m = and m — 2 slow continuum have the 
same frequency 

The plot shown in Fig. [TT] contains the PHOENIX results (crosses) together with the solutions of the four-mode 
coupling scheme represented by Eq. (jUTj) (colored lines), and illustrates the excellent agreement between the two methods. 
As for the cool prominences, the m = 1 Alfven continuum (red) intersects the m = (blue) and m = 2 slow continuum 
(yellow). This may in turn lead to mode coupling, creating a high frequency Am = 1 gap. The analysis of these Am = 1 
gaps was not performed here, but left to future work. We note that the crosses with zero frequency correspond to the 
Eulcrian entropy continuum. 

Fig. ll2l contains the plot for an equilibrium with intermediate gravity g = 0.500. The temperature is again assumed to 
be a flux function. Near the q = 1 surface (s ps 0.927), the results of PHOENIX and the four-mode coupling scheme show 
excellent agreement for the m = 1 Alfven and slow continuum. Owing to the non-zero gravity, a gap has been created 
in the m = 1 Alfven continuum. However, m = and m = 2 slow continua of the coupling scheme slightly mismatch the 
PHOENIX results. This mismatch is because by our neglect of the coupling with m = — 1 and m = 3 slow continua. In 
the plot, these continua are represented by the crosses with non-zero frequency in combination without a solid line. At 
four different radial flux coordinates, s ps 0.82, s w 0.837, s rj 0.973, and s w 0.978 a Am = 1 gap was produced (note 
the color change of the solid lines about these radii). As expected, the four-mode coupling scheme only partly captures 
this kind of coupling. Therefore, about these Am = 1 gaps there is a mismatch between the results of both methods, 
and a separate detailed analysis of the Am = 1 coupling schemes appears to be required. 
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Fig. 10. A part of the continuous MHD spectrum of double layered solar prominence with gravity g = 0.001 computed 
by PHOENIX as a function of the radial flux coordinate s = y/tp. 
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Fig. 11. The real part of the sub-spectrum of the MHD continua of a double layered solar prominence equilibrium with 
small gravity (g = 0.001) as a function of the radial flux coordinate s = are shown for axial wavenumber k = — 1 
and poloidal mode numbers m = {— 2, 4}. Here, the temperature is assumed to be a flux function. The crosses are the 
PHOENIX results, while the colored lines arc the results of the four-mode coupling scheme. The color red, green, blue, 
and yellow correspond to the dominant rj l , £ l5 £ , and C 2 component of the eigenvector. 

The last equilibrium of this subsection that we discuss is the one with strong gravity g = 1.000. The temperature is 
assumed to be a flux function. For this equilibrium, one ob serves the onset of the doub le layered structure in the pressure 
and density as demonstrated in our accompanying paper (jBlokland fc Keppensll2""01lf ). The results for the m = 1 Alfven 
and slow continuum show excellent agreement near the q — 1 surface (s ~ 0.927). Owing to the stronger gravity compared 
to the previous case, the gap in m = 1 Alfven continuum is larger. Near this surface, the previously mentioned mismatch 
of both the m = and m = 2 slow continuum becomes larger because of the neglect of the coupling with both the 
m = — 1 and m = 3 slow continuum. This coupling becomes stronger if the gravity is increased. The four- mode coupling 
scheme is also unable to capture the coupling at the radial flux coordinates s e [0.828, 0.833] and s G [0.975, 0.978]. This 
is again because we neglected the coupling with both the m = — 1 and m = 3 slow continuum. Therefore, close to these 
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Fig. 12. The real part of the sub-spectrum of the MHD continua of a double layered solar prominence equilibrium with 
strong gravity (g = 0.500) as a function of the radial flux coordinate s = \f^j are shown for axial wavenumbcr k = — 1 
and poloidal mode numbers m = {—2, 4}. Here, the temperature is assumed to be a flux function. The crosses are the 
PHOENIX results, while the colored lines are the results of the four-mode coupling scheme. The color red, green, blue, 
and yellow correspond to the dominant rj 1 , f 1; Co^ an d C2 component of the eigenvector. 
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Fig. 13. The real part of the sub-spectrum of the MHD continua of a double layered solar prominence equilibrium with 
strong gravity (g = 1.000) as a function of the radial flux coordinate s = are shown for axial wavenumber k = — 1 
and poloidal mode numbers m = {—2, 4}. Here, the temperature is assumed to be a flux function. The crosses are the 
PHOENIX results, while the colored lines are the results of the four-mode coupling scheme. The color red, green, blue, 
and yellow correspond to the dominant rj 1 , f l5 £ , and £ 2 component of the eigenvector. 



radii there is a large discrepancy between PHOENIX and the coupling scheme in Eq. (|6ip . Furthermore, we note the 
coupling between these two slow continua at the q = 1 surface. This kind of coupling creates a Am = 2 gap that was 
neglected in our theoretical analysis. 

As for the equilibrium class of cool prominences, we observe that when the gravity is increased, the gap in the m = 1 
Alfven continuum increases. This gravity dependence of the gap was investigated in detail by computing the Alfven 
frequency at the 9 = 1 surface for different values of the gravity. The results of this investigation are presented in Fig. 1141 
The results of PHOENIX and the four-mode coupling scheme show excellent agreement with each other. Furthermore, 
we note that the frequency varies almost linearly with the gravity parameter g. This linear dependence is very different 
from the clear non-linear dependence found for cool prominences. This excellent agreement is again only valid for the 
Alfven continuum. As demonstrated in Figs. [T2l and [T3l this does not automatically mean that it also holds for the slow 
continua. The dimcnsionless gravity parameter g — 0.001, g — 0.500, and g = 1.000 correspond to rpog/po = 0.0009, 
rpog/po = 0.4597, and rp g/p = 0.9192, respectively. 
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Fig. 14. The frequency of the Alfven continuum mode rj 1 at the resonant q = 1 surface as a function of the gravity. The 
temperature is assumed to be a flux function. 

7.2. Continua for strong gravity, with axial flow 

We now discuss the last equilibrium of the second class. For this equilibrium, we assume that the temperature is a 
flux function, we adopt an extra strong gravity parameter value (q = 5.000), and we include an axial, shear flow. The 
derivation in our accompanying paper (|Blokland fc Keppensll2"bll[ ) shows that the axial flow must be a flux function. 
The flow itself is described by the profile 

v z (iP)=A i (l-iP), (68) 

where the constant A4 = 0.11. This corresponds to a maximu m Alfven Mach number of 0.105. Using the typical Alfven 
speed of 218 kms -1 , as mentioned in our accompanying paper (jBlokland &: Keppensll20l"i[) . this translates into a filament- 
aligned flow of at most 22.9 km s _1 . Part of the resulting continuous MHD spectrum computed by PHOENIX is shown 
in Fig. 1151 Instead of using 7 poloidal harmonics, we used 15 harmonics to obtain an accurate result (i.e., one that no 
longer changed when adding even more poloidal harmonics). This is a clear indication that the strong gravity causes 
strong coupling between different continua. Owing the extra strong gravity, there is no meaningful comparison between 
the PHOENIX results and those predicted from our theoretical analysis, which is only appropriate if the gravity is 
weak. In the figure, we indicate the non-zero Eulerian entropy continuum (where uj = kv z labeled E) that breaks the 
up-down symmetry about the zero frequency. There is instead now symmetry about this Eulerian entropy continuum. 
Furthermore, PHOENIX finds modes with zero frequency, which correspond to the solution of the generalized eigenvalue 
equation in combination with the perturbed magnetic field being divergence free. If one carefully examines the plot, 
one notices that one continuum branch intersects the Eulerian entropy continuum twice, namely at s 0.808 and at 
s ~ .917. At these locations exactly, the safety factor q is one. Fig. 14 of our accompanying paper (jBlokland fc Keppensl 
|2011|) clearly shows the q = 1 surface locations, as seen in a plot of safety factor q across the vertical direction. Using 
the radial flux coordinate s as a flux surface label where we go from the magnetic axis to the flux rope edge, there will 
be two q = 1 surfaces. We note that the clear coupling effects at other radial locations as well, creating gaps that may 
define preferred frequency ranges for global cigenoscillations. 

8. Conclusions 

We presented the equations describing the complete MHD spectrum. These equations have been specialized for the con- 
tinuous MHD spectrum using techniques frequently used in research on fusion plasma oscillations. From the equations for 
the continuous spectrum, a stability criterion was derived, showing that prominence equilibria for which the temperature 
or the entropy is a flux function have stable continuum modes. However, when the density is a flux function, the continua 
may turn unstable resulting in the convective continuum instability. Whether this instability is related to the intricate 
dynamics witnessed in actual prominence plasmas, or explains why certain filaments disappear abruptly, remains to be 
investigated by follow-up nonlinear computations. We have demonstrated how the onset of CCI modes in flux ropes, which 
carry embedded prominences, will manifest itself as truly flux-surface localized fluctuations. Their growth rate increases 
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Fig. 15. A part of the continuous MHD spectrum of a double layered solar prominence with strong gravity g = 5.000 
computed by PHOENIX as a function of the radial flux coordinate s = \f^. 



when a flux rope expands (the gravity parameter g increases because of both rope widening and the associated mag- 
netic field weakening) or when it becomes more mass-loaded. We note that this route to potentiall y violent prominence 
dyna mics is here related to internal prominence properties. This contrasts with studies (see e.g. MacTaggert fc Hoodl 
l2009f) that emphasize the role of external influences, studies that emphasize the response to changes in the overlying 
field, or the associated pressure buildup in the surroundings by the shearing of its footpoints. We note that the nonlinear 
simulations that may ultimately decide between these competing or possibly cooper ating aspects, will need to go be- 
yond the zero beta, and/or no gravity assumption typically employed in these studies ([Aulanier et al.ll2002t lAmari et al.l 
Il999f) . Moreover, to capture the details of flux-surface localized continuum modes with varying MHD character (slow 
versus Alfvcn) from one flux surface to another, the requirements of the 3D resolution will likely only be achievable by 
using effectively high resolution, grid-adaptive simulations. In additi on, guidance on the ac tual magnetic field structure 
in prominence surroundings will be needed (such as that obtained in lSchmieder et afll2007l) . to distinguish between the 
theoretical freedom in flux profiles. 

The equations for the MHD continuous spectrum have been specialized for small gravity, using an expansion in 
poloidal Fourier harmonics. The resulting equations have been applied to a low frequency Am = gap near a resonant 
q surface. Near this surface, a four-mode coupling scheme has been identified. This scheme couples one Alfvcn mode and 
three slow modes, and can cause gaps (avoided crossings) between continuum branches as these are traced across all 
nested flux surfaces. These gaps appear at low frequency, and form the ideal locations for global discrete eigenoscillations. 
It remains to be shown which kind of coupling occurs at finite frequency, Am = 1 gaps, and to determine whether discrete 
global modes exist within these frequency ranges, under realistic prominence conditions. These global modes, as well as 
their coupling to flux localized surface modes (with resonant damping), need further study, to address whether previously 
established prominence seismology results are consistent with these more realistic prominences. 

The MHD spectral code PHOENIX has been used to accurately compute the continuous MHD spectrum of prominence 
equilibria. Two classes of equilibria have been considered, cool prominences embedded in hot medium and prominences 
with a double layered pressure and density structure for differ ent values of the dimensionless gravity parameter as 
discussed in our accompanying paper (jBlokland &: Keppensll20l"l|) . The results of PHOENIX are compared with analytic 
results of the four- mode coupling scheme, and we have found excellent agreement for small gravity (as expected), while 
for strong gravity, we have identified progressive discrepancies, emphasizing the need for full numerical diagnosis. For 
cool prominences, this discrepancy manifested itself in the slow continuum with the same poloidal mode number as the 
Alfven continuum, while for double layered prominences, differences occurred in the slow continuum with a poloidal mode 
number ±1 with respect to the Alfvenic one. For both classes, when temperature is a flux function, the results demonstrate 
that a gap in the Alfven continuum can appear due to the gravity. For cool prominences, the Alfven frequency associated 
with the gap depends non-linearly on the gravity parameter, while for double layered prominences it scales linearly with 
it. 
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For cool prominences for which the density or the entropy is a flux function, the continuous MHD spectrum has also 
been computed. In the case of the density, the continuous spectrum becomes unstable leading to convective continuum 
instability. Our analysis reveals that close to the q = 1 surface, the CCI mode has a dominant Alfven character, which 
changes away from this surface to a dominant slow type. Furthermore, no gap in the Alfven continuum appears, in 
contrast to the case when the temperature is a flux function. For the entropy case, the continuous spectrum is stable 
and again no gap appears in the Alfven continuum. In the context of prominence configurations, long-lived equilibria 
almost necessitate the temperature flux function dependence, but for active region prominences, the findings about CCI 
possibilities may provide important information about the onset of turbulent dynamics. 

Acknowledgements. This work was carried out within the framework of the European Fusion Programme, supported by the European 
Communities under contract of the Association EURATOM/FOM. Views and opinions expressed herein do not necessarily reflect those 
of the European Commission. RK acknowledges financial support by project GOA/2009/009 (K.U.Leuven). The research leading to these 
results has received funding from the European Commission's Seventh Framework Programme (FP7/2007-2013) under the grant agreement 
SWIFF (project nr 263340, www.swiff.eu). 



References 

Amari, T., Luciani, J., Mikic, Z., & Linker, J. 1999, ApJ, 518, L57 
Arrcgui, I. & Ballester, J. 2011, Space Sci. Rev., in press 
Aulanier, C, De Vore, C, & Antiochos, S. 2002, ApJ, 567, L97 

Belien, A., Botchcv, M., Goedbloed, J., van der Hoist, B., & Keppens, R. 2002, J. Comp. Phys., 182, 91 
Belien, A., Poedts, S., & Goedbloed, J. 1997, A&A, 332, 995 
Berger, T., Shine, R., Slater, C, et al. 2008, ApJ, 676, L89 

Bernstein, I., Frieman, E., Kruskal, M., & Kulsrud, R. 1958, Proc. Roy. Soc. of London, A, 17 

Blokland, J. & Keppens, R. 2011, A&A, (submitted) 

Blokland, J., Keppens, R., & Goedbloed, J. 2007a, A&A, 467, 21 

Blokland, J., van der Hoist, B., Keppens, R., & Goedbloed, J. 2007b, J. Comp. Phys., 226, 509 

Bondeson, A., Iacono, R., & Bhattacharjee, A. 1987, Phys. Fluids, 30, 2167 

Chance, M., Green, J., Grimm, R., & Johnson, J. 1977, Nucl. Fusion, 71, 65 

DeRosa, M., Schrijver, C, Barnes, G., et al. 2009, ApJ, 696, 1780 

Frieman, E. & Rotenberg, M. 1960, Reviews of Modern Physics, 32, 898 

Goedbloed, J. 1971, Physica, 53, 501 

Goedbloed, J. 1975, Phys. Fluids, 18, 1258 

Goedbloed, J. 1997, Transactions of fusion technology, 33, 115 

Goedbloed, J., Belien, A., van der Hoist, B., & Keppens, R. 2004a, Phys. Plasmas, 11, 4332 
Goedbloed, J., Belien, A., van der Hoist, B., & Keppens, R. 2004b, Phys. Plasmas, 11, 28 

Goedbloed, J., Huysmans, G., Holties, H., Kerner, W., & Poedts, S. 1993, Plasma Phys. Control. Fusion, 35, B227 

Goedbloed, J., Keppens, R., & Poedts, S. 2010, Advanced magnetohydrodynamics (Cambridge University Press) 

Goedbloed, J. & Poedts, S. 2004, Principles of magnetohydrodynamics (Cambridge University Press) 

Goossens, M., Erdelyi, R., & Ruderman, M. 2011, Space Sci. Rev., in press 

Hain, K. & Lust, R. 1958, Z. Naturforsch., 13a, 936 

Hameiri, E. 1981, J. Math. Phys., 22, 2080 

Hameiri, E. 1983, Phys. Fluids, 26, 230 

Hellsten, T. & Spies, G. 1979, Phys. Fluids, 22, 743 

Keppens, R., Casse, F., & Goedbloed, J. 2002, ApJ, 569, L121 

Keppens, R., van der Linden, R., & Goossens, M. 1993, Solar Phys., 144, 267 

Labrosse, N., Heinzel, P., Vial, J.-C, et al. 2010, Space Sci. Rev., 151, 243 

Lin, H., Penn, M., & Kuhn, J. 1998, ApJ, 493, 978 

Lin, Y., Engvold, O., Rouppe van der Voort, L., Wiik, J., & Berger, T. 2005, Solar Phys., 226, 239 
Luna, M., Terradas, J., Oliver, R., & Ballester, J. 2010, ApJ, 716, 1371 

Mackay, D., Karpen, J., Ballester, J., Schmieder, B., & Aulanier, G. 2010, Space Sci. Rev., 151, 333 

MacTaggert, D. & Hood, A. 2009, A&A, 508, 445 

Pao, Y. 1975, Nucl. Fusion, 15, 631 

Poedts, S. & Goossens, M. 1991, Solar Phys., 133, 281 

Poedts, S., Hermans, D., & Goossens, M. 1985, A&A, 151, 16 

Poedts, S. & Schwarz, E. 1993, J. Comp. Phys., 105, 165 

Rappaz, J. 1977, Numer. Math., 28, 15 

Schmieder, B., Gunar, S., Heinzel, P., & Anzer, U. 2007, Solar Phys., 241, 53 
Sleijpen, G. & van der Vorst, H. 1996, SIAM J. Matrix Anal. Appl., 17, 401 
Soler, R., Arregui, I., Oliver, R., & Ballester, J. 2010, ApJ, 722, 1778 
Soler, R., Oliver, R., & Ballester, J. 2008, ApJ, 684, 725 
Tandberg-Hanssen, E. 2011, Solar Phys., 269, 237 

Terradas, J., Molowny-Horas, R., Wiehr, E., et al. 2002, A&A, 393, 637 

Trujillo Bucno, J., Landi Degl'Innoccnti, E., Collados, M., Merenda, L., & Manso Sainz, R. 2002, Nature, 415, 403 
van der Hoist, B., Belien, A., & Goedbloed, J. 2000a, Phys. Plasmas, 7, 4208 
van der Hoist, B., Belien, A., & Goedbloed, J. 2000b, Phys. Rev. Lett., 84, 2865 



